Subdivision bifiltration

In topological data analysis, a subdivision bifiltration is a collection of filtered simplicial complexes, typically built upon a set of data points in a metric space, that captures shape and density information about the underlying data set. The subdivision bifiltration relies on a natural filtration of the barycentric subdivision of a simplicial complex by flags of minimum dimension, which encodes density information about the metric space upon which the complex is built. The subdivision bifiltration was first introduced by Donald Sheehy in 2011 as part of his doctoral thesis (later subsumed by a conference paper in 2012) as a discrete model of the multicover bifiltration, a continuous construction whose underlying framework dates back to the 1970s. In particular, Sheehy applied the construction to both the Vietoris-Rips and Čech filtrations, two common objects in the field of topological data analysis. Whereas single parameter filtrations are not robust with respect to outliers in the data, the subdivision-Rips and -Cech bifiltrations satisfy several desirable stability properties.

Definition

Let $T$ be a simplicial complex. Then a nested sequence of simplices $\sigma _{1}\subset \sigma _{2}\subset \cdots \subset \sigma _{k}$ of $T$ is called a flag or chain of $T$. The set of all flags of $T$ comprises an abstract simplicial complex, known as the barycentric subdivision of $T$, denoted by $\operatorname {Bary} (T)$. The barycentric subdivision is naturally identified with a geometric subdivision of $T$, created by starring the geometric realization of $T$ at the barycenter of each simplex.

There is a natural filtration on $\operatorname {Bary} (T)$ by considering for each natural number $k$ the maximal subcomplex of $\operatorname {Bary} (T)$ spanned by vertices of $\operatorname {Bary} (T)$ corresponding to simplices of $T$ of dimension at least $k-1$, which is denoted ${\tilde {\mathcal {S}}}(T)_{k}$. In particular, by this convention, then ${\tilde {\mathcal {S}}}(T)_{1}=\operatorname {Bary} (T)$. Considering the sequence of nested subcomplexes given by varying the parameter $k$, we obtain a filtration on $\operatorname {Bary} (T)$ known as the subdivision filtration. Since the complexes in the subdivision filtration shrink as $k$ increases, we can regard it as a functor ${\tilde {\mathcal {S}}}(-):\mathbb {N} ^{\operatorname {op} }\to \mathbf {Simp}$ from the opposite posetal category $\mathbb {N} ^{\operatorname {op} }$ to the category $\mathbf {Simp}$ of simplicial complexes and simplicial maps.

Let $P$ be a partially ordered set. Given a simplicial filtration $F:P\to \mathbf {Simp}$, regarded as a functor from the posetal category of $P$ to the category $\mathbf {Simp}$, by applying the subdivision filtration object-wise on $F$, we obtain a two-parameter filtration ${\mathcal {S}}(F):\mathbb {N} ^{\operatorname {op} }\times P\to \mathbf {Simp}$, called the subdivision bifiltration.

In particular, when we take $F$ to be the Rips or Čech filtration, we obtain bifiltrations ${\mathcal {S}}\operatorname {Rips} (-)$ and ${\mathcal {S}}\operatorname {{\check {C}}ech} (-)$, respectively.

Properties

The subdivision-Čech bifiltration is weakly equivalent to the multicover bifiltration, implying that they have isomorphic persistent homology. A combinatorial proof of this statement was given in Sheehy's original conference paper, but a more algebraic version was presented in 2017 by Cavanna et al. The ideas from Cavanna's proof were later generalized by Blumberg and Lesnick in a 2022 paper on 2-parameter persistent homology.

By the size of a bifiltration, we mean the number of simplices in the largest complex. The subdivision-Čech bifiltration has exponential size as a function of the number of vertices. This implies that its homology cannot be directly computed in polynomial time. However, for points in Euclidean space, the homology of subdivision-Čech can be computed in polynomial time, up to weak equivalence, via a construction known as the rhomboid bifiltration. As a precursor to the rhomboid bifiltration, Edelsbrunner and Osang presented in 2021 a polyhedral cell complex called the rhomboid tiling, which they used to compute horizontal or vertices slices of the multicover bifiltration up to weak equivalence. This was extended a year later by Corbet et al. to the rhomboid bifiltration, which is weakly equivalent to the multicover bifiltration, but has polynomial size.